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ABSTRACT 

We analyse the influence of rotation on shapes of pulse profiles of fast-rotating (mil- 
lisecond) pulsars. Corotation has two opposing effects: 1) the caustic enhancement of 
the trailing side (TS) by aberration and retardation (AR), which squeezes the emission 
into a narrower phase interval; 2) the weakening of the TS caused by the asymmetry 
of curvature radiation about the dipole axis. Analysis of the radii of curvature of elec- 
tron trajectories in the inertial observer's frame (lOF) enables these two effects to be 
considered together. We demonstrate that for dipolar magnetic field lines on the TS 
there exists a 'caustic phase' beyond which no emission can be observed. This phase 
corresponds to the zero (or minimum) curvature of the lOF trajectories and maxi- 
mum bunching of the emission. The maximum gradient of polarisation angle (PA) in 
the S-shaped PA curve is also associated with the curvature minimum and occurs at 
exactly the same phase. The asymmetry of trajectory curvature with respect to the 
dipole axis affects the curvature emissivity and the efficiency of pair production, sug- 
gesting a minimum at the caustic phase. Emission over a fixed range of altitudes, as 
expected in millisecond pulsars, leads to broad leading profiles and sharp peaks with 
a cutoff phase on the TS. We apply our results to the main pulse of the 5 ms pulsar 
J1012-f5307. 

Key words: pulsars: general ~ pulsars: individual: J1012-I-5307 - Radiation mecha- 
nisms: non-thermal. 



1 INTRODUCTION 

We analyse two specific issues that are important for the ap- 
pearance of pulse-averaged profiles of pulsars: 1) the radial 
extent of a region emitting at a given radio frequency, and 
2) the asymmetry about the magnetic axis of the curvature 
of electrons' trajectories measured in the observer frame. 
The influence of both of these factors on profile shape is 
closely associated with the rotation of the magnetosphere, 
so they should be most noticeable for rapidly rotating ob- 
jects such as the millisecond pulsars (MSPs). Unlike nor- 
mal pulsars, MSPs do not exhibit radius-to-frequency map- 
ping (RFM, see figs. 16 and 17 in Kramer et al. 1999). 
This might possibly be interpreted as near-surface emission 
with negligible radial extent (McConnell et al. 1996). How- 
ever, Dyks, Rudak & Demorest (2009) have shown that at 
least some components in radio pulse profiles are produced 
when our sightline cuts through the curvature emission from 
radially-extended plasma streams. In this stream-cut sce- 
nario, the observed location of a component is determined 
by the frequency-independent moment at which our line of 
sight crosses the plane of a thin stream. When this happens, 
the component becomes simultaneously recorded at all fre- 
quencies that are emitted at the altitude of the cut. There is 



no RFM, because positions of components are fixed by the 
geometry of the cut. 

This interpretation requires that a broad range of ra- 
dio frequencies is emitted by the stream at a fixed alti- 
tude. Such a possibility arises naturally in the stream-cut 
scenario, because density of plasma in a stream decreases 
in the transverse direction from the stream axis. There- 
fore, regardless of the altitude that is considered, a broad 
range of plasma densities (and plasma frequencies) is avail- 
able at different transverse distances from the stream axis. 
It is then possible to have a radially-extended emission re- 
gion (a stream) which provides a range of plasma frequen- 
cies at each altitude (the 'radius-to-frequency' mapping is 
replaced with the 'stream-diameter-to-frequency' mapping). 
Thus, the radially-extended emission at a fixed radio fre- 
quency is possible without violating the foundation of the 
RFM idea, that relates the observed frequency to the local 
plasma frequency at the emission pointy 

The effect of rotation on a radially-extended source can 

^ Needless to say, even in the stream the density must also de- 
crease with altitude, so the RFM is present, but becomes a second- 
order effect. 
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be expected to make intrinsically symmetric pulse profiles 
appear asymmetric with respect to their centre, changing 
the phases of their components, their intensities and their 
widths. In the following section ([2]) we introduce a simple 
(small-angle) analytical model to determine the likely ob- 
servable effects in the profiles of fast-rotating pulsars with 
emitting regions of finite height. These caustic effects are 
well studied numerically and are routinely employed to in- 
terpret high-energy pulse profiles of pulsars (Morini 1983; 
Smith 1986; Romani & Yadigaroglu 1995; Cheng et al. 2000, 
Dyks & Rudak 2003). Here, however, we focus on the an- 
alytical study of the phenomenon, which is feasible for low 
emission altitudes (r < 0.1-Ric) or small distances from the 
dipole axis. This limitation therefore allows possible appli- 
cations to the radio profiles as well as to specific gamma-ray 
peaks, eg. those generated by the trailing caustics of the 
two-pole caustic model. 

The results are used to create artificial but instructive 
pulsar profiles in Section ([3)|. In Section Q we analyse the 
curvature of the particle trajectories and examine the impact 
of considering curvature radiation on model profiles. Then, 
in Section ([6| , we see how far our results can model the main 
pulse profile of J1012-I-5307. Finally, conclusions are drawn 
in Section ((T]). 



2 THE PHASE-DEPENDENT EFFECTS OF 
ROTATION ON PULSAR PROFILES 

2.1 Derivation of caustic coordinates 

Our 'simplest case' model for an averaged pulse profile is 
limited to the plane of the rotational equator ((^ ~ 90°, 
where ( is the viewing angle measured from the rotation 
axis f2), and assumes the orthogonal inclination of dipolar 
magnetic field {a ~ 90°). The B-field in the CF is assumed 
to have the shape of the static dipole. 

The radio-emission region is assumed to have the form 
of a flaring tube/funnel with negligible thickness. It is radi- 
ally extended, but limited to a single value of the footprint 
parameter s and it is symmetric with respect to the {fl, fl)- 
plane (the fiducial plane) . 

In the CF the radiation from the radial position r is 
emitted at the angle 



(1) 



with respect to the dipole axis p, where the footprint pa- 
rameter s = sinS/sinSio, 9 is the angle between r and pi, 
6io ~ {r/RicY^^ is the 6'-coordinate of the last open field 
lines at the same radial distance r. 

The radiation emitted on the leading side (LS) of the 
dipole axis is therefore detected at the phase 



S r „ r 

--S,/ 2 h (pf 

2 V^ic Rlc 



(2) 



where the term 2r/R\c accounts for the leftward shift caused 
by AR effects, ie. the aberration (contributing one r/Ric) 
and the retardation (contributing another r/Ric). The sym- 
bol (j>f denotes the absolute fiducial phase, determined by 
the detection of a photon emitted from the center of the 



NS exactly at the moment when the dipole axis was in the 
Sl-observer plane (see Dyks 2008 for details). 

The radiation on the trailing side (TS) of the dipole 
axis is detected at the phase: 



Ru 



(3) 



The opposing signs of the r-dependent terms in eq. ^ mean 
that on the TS the AR effects tend to cancel the effect of the 
curvature of B-field lines, leading to specific caustic effects. 
In fact, for arbitrarily large r it is easily possible for the TS- 
radiation to arrive/be observed on the leading side of the 
profile (especially if s is small). 

Let us consider electrons moving upward along a mag- 
netic field line with parameter s on the TS of the open re- 
gion. As their altitude increases, their radiation is detected 
at progressively later phases until, at an altitude of Test, an 
increment in height of an electron's trajectory yields an in- 
crease in the linear term in eq. Q which is greater than 
that of the square root term. Emission is then detected at 
the phase ^cst, which is 'the latest-possible detection phase' 
for a magnetic fieldline with a given s, because further in- 
crease in r does not produce any advance in detection phase. 
For still larger altitudes, the radiation becomes detected at 



phases earlier than 



At phases satisfying 



< 



it is 



therefore possible to observe radiation from particles at two 
different emission altitudes on the same fieldline (because 
the quadratic eq. (|3} for r has two real roots). 

For a given s the height rest and phase ^cst can be found 
by using eq. Q to solve the equation: 
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which gives the magnetospheric location of the optimally 
caustically-enhanced emission: 

9 2 



rest 

-Ric ~ 64* ■ 

Then, from eq. 
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(6) 



Thus, for a ~ ~ 90° no tangent-to-B radiation emitted 
within the open fieldline region (s ^ 1) can be detected later 
than 16° after the fiducial phase, regardless of the emission 
altitude. 

Hereafter the 'caustic' coordinates of eq. ([5]) will be de- 
noted by Test and Scst = Scst{r), and the 'caustic enhance- 
ment phase' by ^cst. 



2.2 Caustic phase compression 

To estimate the deforming effects of aberration and retarda- 
tion on a pulsar's profile, let us assume that all radiation is 
emitted between radial distances ri and r2 (where ri < r2), 
and that Ar < Ric- The radial limits will in general refer to 
some selected frequency, so strictly we should write r^,i and 
rv,2, but the index u will be omitted for simplicity. 

From eq. ([5} the LS radiation emitted between ri and 
r2 is therefore spread over the phase range: 



A^L = 0L(r2) - 0L(r-i) 
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The spreading is amplified by the second bracket of eq. ([7]), 
which increases | A</)l| with respect to the case with the aber- 
ration and retardation ignoredlj Similarly, for the equivalent 
fieldline s on the trailing side we use eq. ([3)l to give the net 
change in phase: 



AtpT = (/!>T(r2) 
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In contrast to eq. the terms of eq. Q have opposite 
signs which makes |A0t| < |A0l| whenever r2 > ri and 
s 7^ 0. Because of this, the radiation on the TS 'piles up' 
within a narrower phase range than on the LS and becomes 
intensified. 

If we interpret eqs. ([7)l and ([8} as component widths 
generated by the same fieldline s on the LS and TS respec- 
tively, then we may compare them by expanding eqs. ((Tjl and 
((8| about any emission radius r with r2 > r > ri, giving 
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for the TS. Since both A(/)l and A^t are to first order lin- 
ear in Ar, the profile asymmetry does not depend on Ar 
(the measure of the asymmetry is the ratio of fiuxes on the 
TS and LS: Rp = Ft/Fl = | A(^l|/| A0t|, see eq.IUbelow). 
Thus, assuming the emissivity in the CF is roughly symmet- 
ric with respect to the (f2, /i)-plane, one cannot change the 
relative intensity of the LS and TS by adjusting the emitter's 
radial extent Ar. 

However, this calculation fails in a direct vicinity of the 
caustic coordinates. If the rest corresponding to s lies be- 
tween ri and r2 then eqs. ^ and pO[) do not give the ob- 
served TS component width (because (f>T in eqs. (|10|l is not 
a monotonic function of r when r is close to rest)- In fact 
eqs. (O and (|1U|I understate the component width because 
the phase shifts above and below rest partly overlap and 
cancel one another. 

Suppose, for a given s on the TS, the radiation below 
rest generates a component width (A^ti) which is at least 
partly overlaid by the component width {A(f>T2) generated 
above rest- Then the observed component width A0t will 
be the larger of A(/)ti and A(/)t2 and both must anyway 
extend to the phase (j)cst{s). To see this in eq. (|8]) we set 
Ari = rest ~ ri [i = 1,2) and expand A^t^ about rest to 
the second order to give: 
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^ Throughout this paper the phase intervals that correspond to 
altitude differences are always defined as the upper-altitude phase 
minus lower-altitude phase. Therefore, Aij>L is negative. 



The first term, by definition (eq.[5}, is zero and the remain- 
der has a sign depending on i, giving 



A0T: 



(-1)' Ric (Ar. 
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This expression is independent of s, except indirectly 
through Ar^. Thus the minimum component width is at- 
tained at Ari = Ar2 = ''^"''^ , where the upper and lower 
phase shifts exactly overlap. 

The results of this section are likely to have impor- 
tant implications for understanding MSP profiles because, 
as mentioned earlier, these pulsars are thought to have deep 
emission zones (i.e. large Ar/Ric) and hence caustic effects 
can be expected at a wide range of phases <j)cst , creating nar- 
row, bright components on the TS in contrast to the broad 
dimmer features on the LS. In any pulsar of given period and 
with fixed ri and r2, there will be a definable caustic zone 
in its profile where component narrowing and intensification 
might be expected. Its range is given by 

Ar 

{4>cst)r^ri — {4'cst)r=r2 ~ 2-5—, (13) 

Kle 

Note this result does not depend on the absolute value of 
the emission heights, and will be greatest for pulsars with 
narrow light cylinders. 



2.3 Caustic intensity compression 

The more strongly a given component is compressed by AR 
effects, the larger will be the radio flux observed at the maxi- 
mum of that component. Thus, assuming uniform emissivity, 
the phase-dependent fiux (profile shape) can be estimated 
as: 

where | A(;/>((;/>)| is a width of a component located at pulse 
phase (j>- The quantity F{(l>) represents the flux at the max- 
imum of that component. 

Eq. p4|l will be satisfactory across most of the profile 
(always including the LS), so that we are able to use the 
first-order eqs. (|9]) and (|10|) . However, eq. (|14|l cannot be 
used in a direct vicinity of the caustic coordinates. This is 
because our model assumes an unphysical radiation pattern 
(the delta function along the electron velocity). This leads 
to the infinite value of the phase-resolved flux at = <f)cat- 
The infinity appears because for Ar decreasing down to zero 
around rest, the amount of radiation emitted (proportional 
to Ar), decreases slower than the detection interval A^t oc 
l/(Ar)2 (eq.[T2l). 

Although this result is formally correct, it is unphysical 
because locally the radio waves are emitted into a solid angle 
of finite width Afjjb- Near the caustic coordinates, that is, for 
(ri, r2) = (rest — Ar, rest + Ar) and Ar — > 0, the phase inter- 
val A(I)t becomes negligible in comparison to A(j)t. In such 
a case, the flux observed at = ^cst becomes determined 
by the solid angle of the elementary beam (ocA^^), and by 
Ar. Therefore, the equation [TJ] is not applicable to the caus- 
tic points of the magnetosphere and cannot be used to re- 
liably predict the phase-resolved flux at the caustic phase 

[F{^ = <Pcst)]. 

Such caustic effects have already been numerically ex- 
plored in discussions of the few pulsars found to have high 
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energy profiles. Eq. ([6} is therefore the smaU-angie equiv- 
alent of the phase at which the leading peak of the two- 
pole caustic or slot-gap model occurs (Dyks & Rudak 2003; 
Muslimov & Harding 2004). This leading gamma-ray peak 
is emitted from low altitudes on the trailing side of the polar 
tube and lags the main radio pulse by ~ O.IP, if the clos- 
est pole is viewed at a small impact angle (C — ce)- In this 
case the phase location of this gamma-ray peak is given by 
eq. (|A4p of the Appendix |X] for arbitrary dipole inclinations 



3 THE GENERATION OF ASYMMETRIC 
PROFILES. 

To illustrate the observable effects of AR on a millisecond 
pulsar we generate artificial profiles for a pulsar assumed 
to have a period of 5.25 mijj and to radiate its emission 
between heights of ri — 10^' cm (approximately the neutron 
star surface) and r2 — 2 ■ 10^ cm. Such a pulsar would have 



''1 



= 0.04 and ^ = 0.08. 



Ric = 25 • 10"" cm so that 

In Fig. [T] an intrinsically symmetric 6-component (i.e. 
3-cone) emission region is assumed for 2 different sets of s 
parameters. The first (Fig. [1^) includes an outermost cone 
on s = 1 and has a total width of 55° {— 3y/j^, the open- 
ing angle of the last closed fieldline at the top of the emis- 
sion zone). From eqs. ^ and (|8]) this gives Acjf>L — 15° 
and A(j>T — 3°, i.e. the leadingmost component is 5 times 
broader than the trailingmost. The second (Fig.lTJj) has nar- 
rower cones and the merging of components is more compre- 
hensive on the LS. In both cases, the radial extent causes 
the LS profiles to assume a 'boxy' shape so that even with 
widely spaced cones the stretched components overlap and 
false 'components' can be created (marked 'O' in Fig. [l^). 
It can also be seen that the stretched components do not 
always fully merge with each other, so that a single notch 
can appear (marked 'N' in Fig. [l^). Single notches in boxy 
profiles were described by Cordes (1975) for B1919-(-21. Sim- 
ulated polarisation profiles (not shown) demonstrate that 
the overlap does not lead to any depolarisation, because the 
overlapping components have similar polarisation angle (and 
our code assumes they have the same intrinsic polarization 
degree) . 

By contrast, the TS components are clearly narrowed 
and become separated. From eq.([5]) it follows that strong 
caustic effects will be observed on fieldlines with s between 
0.53 and 0.75 at phases between 4.5° and 9°. This is precisely 
the phase range within which the strongest peaks occur (on 
s = 0.6 in case (a) and s = 0.65 in case (b)). Note that the 
height of the strongest components is limited by the phase 
resolution of our calculation (360 phase bins per period). 

Note also that, although the total width of the profile 
remains unchanged from its CF value (3^-^ — 55°), any 
central emission on s = (on a straight fieldline in the CF) 
would be spread over 6° between 4.5° and 9° preceding the 
fiducial phase. 

In general, for pulsars with trailing components emitted 
near rest the flux at the trailing edge of a profile is then likely 



This typical MSP period will be used throughout this paper to 
enable a comparison with the profile of J1012-I-5307. 
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Figure 1. Modelled pulse profiles for the emission region extend- 
ing radially between 10® and 2 • 10® cm in a pulsar with P = 5.25 
ms and a = f = 90°. This implies a caustic zone in the pro- 
file ranging from 4.5° to 9°. In a the radiation is emitted from 
s = 0.2, 0.6, and 1.0. In b s = 0.25, 0.45, and 0.65. The profiles 
formed without AR effects are shown as dotted lines. In both cases 
the expected enhancement and narrowing of components occurs 
in the caustic zone of the TS. The result is for a fictitious emis- 
sion mechanism with fixed radiative power per each centimeter of 
electron trajectory in the lOF. Zero on the horizontal axis corre- 
sponds to the absolute fiducial phase and the phase resolution is 
1°. 



to drop more steeply than on the LS. When a ~ 90° and 
s ~ 1 on the trailing edge of a pulse, the caustic phase 
coincides with the trailing edge of a profile if the emission 
occurs at rcst/(10® cm) = 1, 2, 3.3, 6.7, 67, 670 for P = 1.5, 
3, 5, 10, 100, and 1000 ms, respectively. Assuming that the 
radio emission altitude does not exceed a few tens of 7?ns we 
expect steep trailing flanks for profiles of interpulsars with 
P < 100 ms. 

Sharp features on the trailing side of profiles indeed 
seem to be a frequent feature among MSPs (J0030-I-0451, 
Lommen et al. 2000; J0613-0200, J1022-I-1001, B1855-h09, 
B1937+21, Kramer et al. 1999; Stairs et al. 1999). In prac- 
tice it remains difficult to disentangle the components of 
an MSP and hence to make estimates of ri and r2. This is 
mainly because the true position of the fiducial phase within 
the pulse is unknown and the pulsar's angle of inclination, 
which affects our calculation through the generalised forms 
of eqs. ((2} and (|3]) (see Appendix A), is not easy to estimate. 
However, as we show below, other consequences of AR may 
be utilised to give further clues. 

3.1 Position angle swing and caustic enhancement 

A simple physical description of the caustic pile-up of ra- 
diation is that electrons at rest for a short while move rec- 
tilinearly towards the observer in lOF. At r < rest their 
lOF trajectory is bent in the backward/trailing direction. 
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whereas at r > r-cst it is bent forward because of the domi- 
nation of the corotation. It is the inflection point of the lOF 
trajectory that is located at (rcst,Scst) [or (rest, finest)]- Thus, 
the region of the caustic pile-up (eq.[5} is simultaneously the 
region of the zero-curvature of electron trajectory in lOF, 
as localized by eq. (3) in Dyks (2008). By noting that 



(15) 



" (r/iiic)i/2' 
one can rewrite eq. ((Sjl in the forms: 



' est 
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^cst ■ 



(16) 



For a laterally-extended emitter the zero (or minimum) cur- 
vature corresponds to the steepest gradient of the polar- 
ization angle swing. The locations of sharp caustic peaks 
in pulse profiles are then predicted (for a radially-extended 
emitter) to occur at the same phase, ie. ^cst = <?!>pa. Indeed, 
eq. (O can be written as: 

■ ' -'^"+0/^0PA, (17) 



(cf. eq. 16 in Dyks 2008). 

The fact that (/!>cst = (/>pa does not always mean that 
caustic enhancement and the PA inflection point will be 
observed at the same phase in a given object. This is be- 
cause the relativistically-delayed S-swing is expected for a 
laterally-extended emitter, whereas the caustic enhancement 
is for the radially-extended emission, in which case the PA 
curve does not have to assume the S-shape. Some examples 
of the PA curves expected for the extreme case of radial ex- 
tent are shown in Dyks, Harding & Rudak (2004), eg. their 
fig. 5 is for As — and fig. 7 for As ~ 0.05. Rapid changes of 
the PA are expected at the caustic peaks and a rather com- 
plicated PA curve, cf. the observed optical polarisation data 
of the Crab pulsar, Slowikowska et al. (2009). Moreover, for 
the specific case of curvature radiation the emissivity van- 
ishes for zero lOF curvature, which can produce a minimum 
of flux at 0cst in the averaged pulse profile (see fig. 3 in 
Blaskiewicz et al. 1991, and the discussion in our Sect. 14.3]) . 
We anticipate the case of i^cst = 0pa to be observable for a 
slab-shaped emitter with the non-negligible lateral and ra- 
dial extent. 



4 CAUSTICS AND EMISSIVITY 

4.1 Radius of curvature of electron trajectories in 
the observer's frame of reference 

The corotation increases the curvature of electron trajec- 
tories on the LS of the MP, whereas it decreases the curva- 
ture on the TS (Ahmadi and Gangadhara 2002; Dyks 2008). 
The question of how this affects radio emissivity is a com- 
plicated one, because it depends on details of the amplifi- 
cation process. One one hand, one can qualitatively argue 
that the smaller curvature (larger curvature radius) makes 
the amplification more efficient, because the electrons fol- 
lowing less-curved B-field interact with narrow beams of ra- 
diation for a longer period of time. On the other hand, for 
the coherent curvature emission from bunches of n charges, 
the radio emissivity is proportional to the noncoherent cur- 
vature emissivity: Fcoh = n^Fcr, where Fcr oc p~^^^^ is the 



noncoherent curvature emissivity of a single electron in the 
low- frequency limit (Buschauer & Benford 1980), and piof is 
the radius of curvature of electron trajectory in the inertial 
observer frame (lOF). This holds at least as long as n is 
independent of piof. 

For the latter mechanism our model for the pulse profile 
shape becomes: 



F((^) oc 



-2/3 
Piof 



Ar 



\A<t>i^)\ 



(18) 



To use it efficiently, we need a convenient method for esti- 
mating Piof. The value of piof can be expressed as: 



Piof 



l^iof 



(19) 



where Siof is the electron acceleration in the lOF. In eq. ([19} 
the electron velocity has been approximated by the speed of 
light c, and any acceleration parallel to the electron trajec- 
tory is thereby neglected. 

The lOF acceleration can be expressed as a sum of two 
vectors: 

Siof — Sb + an ^ — pB + 2r2c Crf,, (20) 
Pb 

where Sb is the CF-acceleration caused by the curvature 
of magnetic lines and an is the acceleration generated by 
the corotation (see Dyks 2008 for a simple derivation; note 
that we still keep the discussion limited to the orthogonal 
case of a = ~ 90°). The symbol ps denotes the radius of 
curvature of B-field lines in the CF: 



(21) 



pB (with a hat) is a unit vector along that radius of curva- 
ture, and 60 is a unit vector in the direction of the rotational 
azimuth. In the small-angle approximation, on the leading 
side in the plane of the rotational equator we have Pb — e^. 
On the trailing side p^ ~ — e^, so that eq. (|20|l becomes 
one-dimensional. 

Thus, for these specific simple cases we correspondingly 
have: 



Piof.L 




Ru 



Ric 
Pb 



(leading side), 



Piof.T — Rh 



3 s 
4 



-1 




= -Ric 


^-2 




Pb 



(22) 



(23) 



(trailing side). 

Interestingly, the only difference between piof and the 
CF-curvature radius of magnetic field lines (eq. 1211) is the 
added number '±2' in eqs. (|22p and (|23p . Not surprisingly, 
apart from the sign issue eq. (|23p is the inverse of dcp-r/dr 
(eq. 131 . This is because the very definition of the curvature 
radius is piof = dl/dd, where dl is an increment of electron 
trajectory and d5 is the angular increment corresponding to 
dl. Because of the small-angle approximation we have dl ~ 
dr, and because of the near-orthogonality we have dS ~ d(j>, 
ie. a change of the direction tangent to electron trajectory 
in lOF is equal to the interval of the pulse phase. Thus, 
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Figure 2. lOF trajectories (thiclc solid lines) of electrons moving in the equatorial plane of orthogonal pulsar (a = 90°). The curves are 
separated by As = 0.1 and refer to electrons that left the star surface when the dipole axis was vertical in the figure (i.e. the fiducial 
instant). lOF trajectories for electrons moving along the dipole axis are marked with dots (merging into very thick solid lines) and 
coincide with the analytical solution (dashed circles of radius centered at the crosses on both sides of the star; the centres are 

located at x/Ri^ = ±0.5 =F (0.5i?Ns/^lc)^ E'nd v/Rlc = TO.S/Jns/^Ic)- Electron trajectories in the CF (B-ficld lines) are shown as dotted 
lines, except for the last open ones, which are thin solid. The outermost circle is the light cylinder. The figure assumes rigid corotation 
of a static-shape dipole out to 0.95iJic. The trajectories shown extend from the leading edge to the trailing edge of the open region and 
have s = 1, 0.9, 0.8, 0.1, 0.0, 0.1, 0.9, 1.0. 



Piof = dl/d6 ~ {d(l>/dr)-^. For a ~ 90° the ratio of the lOF- 
curvature radii on the leading and traihng side can thus 
be conveniently estimated with the help of eqs. (|22|l and 
H23p (the model is extended to the nonorthogonal dipole 
inclinations in Appendix IX)) . 

For MSPs with P ~ a few ms, the value of piof on the 
trailing side of the polar cap surface can be one order of 
magnitude larger than the equivalent value on the leading 
side. For P = 5 ms and s = 1 we have piof.L ~ 5 • 10^ cm, 
Piof.T ~ 2.5 • 10^ cm, whereas the dipolar ps = 8 • 10® cm. 
The asymmetry of piof about the dipole axis is illustrated 
in Fig. [2] where particle trajectories in the lOF are shown. 
It exhibits trajectories from the leading edge to the trailing 
edge of the open region [there are 21 trajectories that start 
from each polar cap with footprint parameter s that goes 



from 1 on the LS through to 1 on the TS in steps of 0.1]. 
The gradual increase of curvature radii towards the TS is 
visible. The dipole axis is marked with merged dots (dark 
line) . 

In fact, at fixed r the curvature is completely symmetric 
about the caustic trajectory, that is, about Scst(''). This can 
be seen by using eqs. (|22|) and ([23} to transform the basic 
equations ^ and ((3]), giving 



2r 



Piof 



(24) 



Then, if cones are formed with the same lOF curvature 
they will be symmetric about the caustic phase at any given 
height r. If a finite range of r is considered the cones will be- 
come "blurred" - but nevertheless this cannot lead to the TS 
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peaks being brought closer together. Fig. |3] shows curvature 
as a function of altitude for a range of field line parameters 
(and compares this with our approximate analytical formu- 
lae for piof and ps.) This makes clear that our small angle 
approximation holds well in the region we are interested in. 

These results imply that with the knowledge of the cur- 
vature radius in the lOF we can transfer the entire discus- 
sion of components' widths and intensities into the lOF, 
without the need to refer to the effects of aberration and re- 
tardation. Specifically, eqs. ((TJ and ((Sjl are equivalent to the 
natural assumption that the component's width is inversely 
proportional to the radius of curvature in lOF: 

Ar 

A<^~— . (25) 

Piof 

This can be immediately proved by comparison of the ap- 
proximate eqs. ((9} and (I10|l with the formulae for piof 
fegs. 1221 and I23 p . The caustic effects of AR can thus be seen 
as simply the consequence of an increase of the curvature 
radius in lOF. 

Our flux model of eq. (|18p can then be written 



F{(t)) oc 



Ar/pio 



= Piof 



(26) 



For the non-coherent, low-frequency curvature radiation, the 
exponent a is equal to —2/3, and the equation (|26|) just 
reduces to F{(j>) oc Pj'^f , which suggests a weak amplification 
on the TS. However, as we shall see in the next subsection, 
more physical considerations are likely to reverse this effect. 

As a consequence, the TS of the profile is subject to two 
potentially competing effects, both dependent on trajectory 
curvature and hence both having maximum effect around 
the caustic zone of the profile. On the one hand, the AR 
bunches the emission of a single fieldline into a narrow phase. 
On the other, if the coherent radio emissivity decreases with 
increasing piof (as the model of eq. [18] assumes), this may 
counteract the caustic amplification. 
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Figure 3. lOF radius of curvature calculated for selected mag- 
netic field lines of a pulsar with P = 5.25 ms and a = 90°. Solid 
lines present exact numerical solution and dashed lines are for 
eqs. I I22I I and II23I I. The numbers give the footprint parameter s. 
Dotted lines (the lowest one for s = 1) present the dipolar ps- 
The horizontal marked with '0.0' is for the dipole axis trajectory 
(piof = and separates the lower curves for the LS from the 

upper lines for TS. 



4.2 Giant pulses at the caustic singularity? 

It is also interesting to consider the implications of caustic 
enhancement for single pulses. The quasi-rectilinear motion 
of electrons in the caustic regions of magnetosphere makes 
favourable conditions for coherent amplification, because the 
interaction between the emitted radiation and particles can 
last for a longer period of time than in the case of a strongly- 
curved trajectory. We speculate that this may give rise to 'gi- 
ant' pulses, whose phase location is determined by eq.[6](see 
Appendix A for the non-orthogonal version). At the same 
time, however, the reduced curvature decreases the emissiv- 
ity and frequency of the curvature radiation. To overbalance 
this effect, the efficiency of the amplification process would 
have to grow with decreasing curvature sufficiently quickly. 
The energy of electrons has also to be larger for the curva- 
ture radiation to extend up to the observed band. 



4.3 Rotation and the profile flux balance 

Eq. (|26|l enables us to consider the profile asymmetry in 
more realistic cases, differing from the basic cases of Fig.l 
where Ar was held constant across the profiles and only AR 
effects were taken into account. Let us assume an emission 



region with the quantity Ar symmetric with respect to the 
main meridian so that the only asymmetric quantity left is 
Piof. We also allow for a non- negligible range of the footprint 
parameter As associated with a given component because 
it is easy to model numerically and can be expected in re- 
alistic cases. The profile asymmetry is also determined by 
the value of the exponent a that describes how the intrinsic 
coherent emissivity depends on piof. In the discussion below 
we assume a = —2/3 for definiteness [non-coherent curva- 
ture radiation (CR) in the low-frequency limit]. 

The effect of non-zero As can be modelled by modifying 
eq. (|26|l so that 

pfofAr 



F(<^) oc 



Aj^/Piof + iy^As' 



(27) 



In the case of negligible Ar (ie. for Ar/piof ^ 
I As in eq. I27|l we have F oc p^^^^, ie. the compo- 
nents on TS of the profile are dimmer than those on the LS. 
This can be seen in Fig. [4^ which is calculated for Ar = 
(ri,r2) = (0.04, 0.05) 7?io and As = (s^in, s^ax) = (0.9,1). 
Fig. 3 in Blaskiewicz et al. (1991) presents a similar case for 
a region completely filling in the open region [As = (0, 1)]. 

In Fig. |4]d the upper radial distance r2 is increased to 
0.12i?ic. This makes the caustic term Ar/piof in eq. (|26|l 
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Figure 4. Simulated pulse profiles that show how the rotational 
asymmetry is affected by changes of Ar = (ri,r2) and As = 
(smini Smax)- In all figures the intrinsic emissivity is proportional 
to Pj^f ''^ (non-coherent CR). In a) As = (0.9, 1), and Ar/Ri^ = 
(0.04,0.05). In b) r2 is increased to 0.f2Ric. In c) r2 = 0.12i?,ic 
and Smin = 0.5. 



more important so that we have F oc p^J^^^ piof and the TS 
component becomes brighter [the component's 'brightness' 
is understood as the flux at the maximum of the compo- 
nent, not the phase-integrated flux]. For As — (0.5, 1) and 
r2 = 0.12i?ic (Fig. Ht), the decline of the intrinsic emissivity 
{Fcoh (X Pi^f again dominates over the caustic amplifica- 
tion and the LS component becomes brighter. In addition to 
the eq. (|27|l this change of Rf can be understood as follows. 
The peak flux on the LS is increased because the emission 
from adjacent s overlaps in phase. On the TS the caustic 
peaks associated with different s do not occur at the same 
phase. The larger value of As makes the components on the 
TS wider, but their peak flux is not enhanced considerably. 
The influence of Ar and As on the pulse shape can therefore 
be summarized in the following two points: 1) the increase 
of Ar makes the LS components wider and the TS compo- 
nents stronger (larger peak flux. Fig. |31d) ; 2) the increase of 
As works in just the opposite way; it makes the LS compo- 
nents stronger and the TS components wider (Fig. [l];). 
When the emissivity is independent of piof (a — 0) it is 
the TS that is always brighter, with the flux ratio deter- 
mined by the rel ative values of Ar and the intrinsic width 
A,^intr = (3/2)^r/i?icAs. 

We therefore conclude that the caustic amplification (oc 
Piof) makes the TS brighter only if the intrinsic emissivity 
does not decrease with increasing piof too fast. For the low- 
frequency CR (a = —2/3) the caustic amplification can still 
dominate over the intrinsic fiux decline (as in our Fig. IJJj) 
but this requires special conditions to be fulfilled (small As, 
large Ar). 

For large Ar and several emission cones, the above con- 



clusions are furthermore complicated by the possible overlap 
of different components on the LS, which can additionally 
increase the peak flux on the LS. In the simulated profile of 
Fig. [TJj the two maxima at = —18° and —14° are created 
by an overlap of emission pattern from three cones. Such an 
inter-component overlapping is not included in eqs. (|26|l and 
(EZl). 

We then find that it is possible (and not unlikely) for 
the LS to become brighter purely because of the intrinsic de- 
pendence of emission mechanism on piof. For a = —2/3 the 
TS can be enhanced only under special circumstances (small 
As, not overlapping components). Narrow and strong trail- 
ing components seem to be a frequent phenomenon among 
MSPs (Kramer et al. 1998; 1999). Moreover, at least some 
exceptions (with apparently brighter LS) seem to be created 
by the scattering phenomenon (see eg. fig. 8 in Kramer et 
al. 1999). This suggests that —2/3 < a < 0, provided we are 
not biased by unrecognized absorption effects. The results 
of this section also apply for the models of high-energy emis- 
sion (even for the polar cap model if applied for pulsars with 
millisecond periods). In the outer gap and slot gap models, 
the peaks in the gamma-ray profiles are generated mostly 
by the caustic effects. They are then formed in the region 
where the lOF radius of curvature is greatly increased in 
comparison to the CF value. It is then absolutely necessary 
to use the lOF curvature radius piof to obtain reliable pulse 
profiles in the numerical simulations of the CR. 



5 PULSARS WITH PHASE-SHIFTED CONES 

Gangadhara & Gupta (2001, hereafter GGOl), noted that 
profiles of selected pulsars have outer cones more shifted to- 
wards early phase than inner cones. Krzeszowski et al. (2009) 
find at least a tendency towards this effect in most pulsars 
in their sample, but they also find some counterexamples 
(e.g. B1831— 04). The forward shift of outer cones was inter- 
preted by GGOl in terms of a lower emission altitude for an 
inner cone. The detection of such a low-altitude inner cone 
is delayed by a decrease in the aberration angle and the in- 
creased propagation time. This interpretation may be called 
the 'cone-altitude effect'. It implies a structured emission re- 
gion with two rings of different size located at two different 
altitudes (and with a core emission region near the surface) . 

Because such a structure has never received a physical 
explanation it is worth considering whether an alternative 
scenario exists in which the cones can be actually located at 
the same altitude, but intrinsically decentered in the CF to 
produce the forward shift of the cones. A possible natural ar- 
gument for this is the shift of the pair-production symmetry 
point towards the TS. 

As shown in Fig. El for P = 5.25 ms and r = 1.3 ■ 10** 
cm, the zero curvature of electron trajectories in the lOF is 
expected at Scst = 0.61 on the TS of the open region, ie. at a 
phase (j) ~ 0MP1 + 0.8Wmp, where Wmp is the profile width 
and ^MPi is the location of the first (leadingmost) compo- 
nent y Conversely, any non-curvature photons emitted by 
the electrons (moving along the _B-field line with s = 0.61) 

* In the case of a rotationally-distorted vacuum dipole the phase 
is equal to 0cst(fiNs) = </'MPi +0.68Wmp, because of the bound- 
ary of the open region is shifted toward the TS. 



© 2002 RAS, MNRAS OOO.nTOI 



The rotational asymmetry of pulsar profiles 9 




-1.0 -0.5 0.0 0.5 1.0 



Figure 5. Location of the zero-curvature point (bullet marked 
'ZC') at the polar cap of a pulsar with a = 90°, P = 5.25 ms, 
Rns = 1-3 ■ 10® cm. The curvature of electron trajectories in lOF 
(thereby the curvature radiation) vanishes at the point where s = 
0.61. Curved CF trajectories of non-curvature photons that are 
emitted there coincide with the local B-field line in the CF. Thick 
solid and dotted lines mark the polar cap edge for the retarded 
and static-shape dipoles. Thin solid circle marks the standard 
circular polar cap (formally for a = 0°). The decentered circles 
of dots show possible locations of emission cones affected by the 
ZC point. 



propagate in the CF along a curved path locally coincident 
with this B-field line and are least likely to participate in 
any pair-creation cascade (see fig. 7 in Dyks & Rudak 2002) . 
The point of zero lOF curvature is therefore also the loca- 
tion of vanishing (or minimum) pair production. The density 
distribution of the created e±-pairs is then expected to be 
symmetrical around the ZC point located clearly on the TS 
of the open region. The same symmetry (centered at the 
zero-curvature coordinates) can be imposed on the acceler- 
ating electric field if it is strongly screened by the pairs. 

On the other hand, the general geometry of the open 
region implies the dipole-centered symmetry. If the outer 
boundary s = 1 has indeed the same electric potential ev- 
erywhere, the (unscreened) electric field should roughly fol- 
low the symmetry of the open dipolar fieldlines (which are 
centered on the dipole axis, not at the lowest-curvature lOF- 
trajectory). It is therefore conceivable that any conal emis- 
sion structures within the open region will be affected by 
both symmetries and may assume the 'intermediate-case' 
geometry shown in Fig. [S] This effect can mimic the cone- 
altitude effect even if all the cones originate at the same 
altitude. 
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Figure 6. The main pulse of J1012+5307. Radial extent of the ra- 
dio emission implies different widths of components on the leading 
and trailing side of the MP. The two pairs of vertical lines present 
the components' widths expected for low-altitude emission (be- 
tween ri = 10'^ and r2 = 2.3 ■ 10"^ cm) from the last open field 
lines (s = 1). 

6 APPLICATION TO THE MAIN PULSE OF 
J1012+5307 

6.1 Component widths 

Fig. [6] presents the main pulse (MP) of a 5.26 ms pulsar 
J1012-I-5307 observed at 0.82 GHz (top) and 1.4 MHz (bot- 
tom) with the Green Bank Telescope. The bandwidths were 
64 MHz at both frequencies and the total integration time of 
~15 hours spanned the period between July 2004 and March 
2007. This pulsar was chosen because it exhibits an inter- 
pulse (not shown) , strongly suggesting a near-orthogonal ori- 
entation of rotational and dipole axis. The main pulse profile 
exhibits asymmetry which looks consistent with the implica- 
tions of eqs. (O and ((Sjl : the components on the leading side 
are very broad and strongly merged with each other, whereas 
the components on the trailing side, despite being located 
much closer to each other, can still be identified through the 
noticeable minima between them. The trailing components 
must then be much narrower than the leading ones. Because 
of the blending the exact width of the components cannot 
be determined. However, because the components are never- 
theless identifiable, we assume that their widths are similar 
to the components' separation. 

For the leadingmost component of the MP we use the 
width of A(t>L — —15° which is marked by the leftmost pair 
of vertical solid lines in Fig. [6] The width can be used to es- 
timate the corresponding range of radial distance Ar. From 
eq. ([7} we have: 




For ri = lO'' cm and s = 1 we get r2 — 2.3 • 10® and 
Ar = r2 — ri ~ 1.3 • 10® cm. Thus, the width of components 
on the leading side of the MP implies the radial extent of 
the order of the NS radius _Rns ■ 
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Assuming approximate symmetry of the outer cone with 
respect to the {Q, /2)-plane in the CF, one can use the de- 
rived value of Ar in eq. (HJ to estimate the width of the 
components on the TS of the MP and thereby to test the 
reasoning. For ri — 10^ cm and r2 = 2.3- 10^ cm eq. Q gives 
A(j)T — 2.9°, which is marked in Fig. [§] by the pair of the 
vertical solid lines on the right-hand side. The width-scale of 
the leading and trailing components is then consistent with 
the radio emitter's radial extent of the order of ~ 10® cm. 

6.2 Component heights 

However, the interpretation shown in Fig. 6 cannot be made 
consistent with the symmetry properties of piof. According 
to eq. (|24p . the lOF radius of curvature is symmetrical with 
respect to the caustic phase. Thus, both the caustic effects 
and piof-driven emissivity should be symmetrical with re- 
spect to 0cst. For r < 0.14_Ric the caustic phase always pre- 
cedes the trailingmost component (if the latter is assumed to 
have s = 1). The height of components, ie. the flux F{(j}), is 
then expected to fall off symmetrically on both sides of some 
phase located within the pulse window. In the main pulse of 
J1012-I-5307 the flux increases monotonically towards its ex- 
treme trailing edge which can be interpreted in two ways: 1) 
the profile is shaped by physical effects independent of piof; 
2) the emission region is not centered at the dipole axis, but 
is located fully on the LS of the caustic coordinates. 

In the second case, (jicst needs to be located somewhere 
on the right-hand side of the trailingmost component. Since 
no observed features can be noticed there, and the available 
model parameters allow for some fit ambiguity, in Fig. [7] 
we present just a sample profile calculated for (j>cst located 
at the phase of the trailingmost component. The modelled 
profile (bottom solid line) has been obtained for ri — 10^, 
r2 = 2 - 10*^ cm. As = 0.1, q = C = 90° and for uniform 
emissivity per centimeter of lOF trajectory. The component 
on the extreme right-hand side is located right at the caustic 
phase (fleet — 7° (for rav = 1.5 • 10® cm). Its width and 
height are therefore mostly determined by the magnitude 
of As. The flux on the LS is dominated by the effects of 
Piof. In these two limits the flux is well approximated with 
the analytical model (|27p (dashed line, symmetric around 
<^cst). In the intermediate region [phase range (—5°, 5°)] the 
agreement is worse. The values of s parameter have been 
selected for each component to match the observed pulse 
phase and vary from s = 1 for the component on the leading 
edge, through s = near the central minimum, to the caustic 
value of s = 0.6 for the trailingmost component. 

The numerical proflle exhibits the general observed 
properties, ie. the monotonic increase of flux towards the 
TS and the narrower components on the TS. The LS/TS 
flux ratio could be made more consistent with data by plac- 
ing (pest far on the traihng side. However, such fit would not 
be unique, because the ratio can also be adjusted by chang- 
ing As and Ar. The difficulty in establishing 0cst makes 
any multidimensional fitting procedure somewhat ambigu- 
ous and discourages us from forcing the model to closely 
match the data. The LS/TS asymmetry can also be prob- 
ably decreased by increasing the impact angle (3 — — a. 
Moreover, the observed profiles seem to be affected by the 
overlapping of neighbouring components. This cannot be 
easily modelled with our analytical formulae nor with the 
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Figure 7. Numerically-calculated pulse shape (bottom curve) 
drawn below the J1012+5307 data. This uniform-emissivity 
model has ri = 10^ cm, r2 = 2 ■ 10^ cm, As = 0.1, and 
a = f = 90°. Dashed curve presents eq. II27I I with a = and 
assumes that (pest = 7° . 

numerical code, which creates pronounced spurious compo- 
nents caused by the sharp boundaries of the emission region 
at ri and r2. To avoid them, the code would have to be 
supplied with additional parameters that control the radial 
profile of emissivity. 

The CR is physically better justified than the uniform 
emissivity model shown in Fig. [71 and it has been shown to 
produce enigmatic bifurcated components in the profile of 
J1012-I-5307 (Dyks, Rudak & Demorest 2009). In the case 
of the CR, a minimum of flux is expected at 0cst . This sug- 
gests that the caustic phase be associated with the central 
minimum of the MP. We do not model this case, however, 
because the observed MP is not symmetric about this min- 
imum (whereas the modelled profile would have been sym- 
metric) . 

An alternative way to successfully model the profile is 
the cone-altitude approach of GGOl, which assumes a set of 
the dipole-axis-centered rings of emission in the CF. Assum- 
ing that the near-surface emission within the main meridian 
is detected at the central minimum, the observed profiles 
can be reproduced if the rings are located at s = 0.53, 0.7, 
and 0.84, and the average emission radius is equal to 1, 1.7, 
and 2.5 ■ 10® cm, respectively. 

Apart from J1012-f 5307, a similar systematic increase 
of flux towards the TS is exhibited by the central structure 
in the 438 MHz profile of J0437-4715 (Navarro et al. 1997). 



7 CONCLUSIONS 

We find that the shape of a profile is affected by two oppos- 
ing effects of rotation: 1) the caustic enhancement of the TS, 
and 2) the increase of piof on the TS which decreases the cen- 
tripetal acceleration and can thereby weaken the emissivity. 
The first can be viewed as a direct result of the straightening 
of an electron trajectory in the lOF, because the radiation 
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emitted within Ar is received within ~ Ar/piof. The ra- 
diation that constitutes the caustic peaks is emitted from the 
inflection points of lOF trajectories (or from regions of min- 
imum lOF curvature in a non-orthogonal case). Formally, 
the caustic peaks are expected to occur at the phase of the 
center of the relativistically-shifted S-curve of the polarisa- 
tion angle. The caustic phase, if calculated for the maximum 
s = 1, becomes the 'maximum possible detection' phase for 
tangent-to-B radiation within the open region. Fast-rotating 
and highly-inclined pulsars are then expected to have pulse 
profiles with sharply defined (distinctively marked) trailing 
sides. On the LS, AR causes radiation from the outflow- 
ing electrons to be detected at a phase that monotonically 
increases with the emission altitude. The LS of the profile 
is thus expected to grow only gradually, creating a leading 
wing (the extent of which in practice is likely determined by 
the cooling timescale of the electrons) . 

Whether the leading or trailing side of the profile is the 
more enhanced by rotation is likely to be determined by the 
intrinsic dependence of the coherent emission mechanism on 
Piof. In the case of the CR the emissivity is entirely driven 
by the macroscopic centripetal acceleration and we may ex- 
pect a minimum in the zero curvature region. Otherwise, 
a caustic peak will dominate the profile (cf. Fig. [T]3 with 
UJ:). Statistical analysis of millisecond profiles, could there- 
fore provide limits on the possible value of the exponent a 
in the relation Fcoh oc /9?,f. Inspection of published profiles 
suggests that a < (as opposed to a = —2/3), although the 
number of good-quality millisecond profiles is very limited, 
and the interpretation of profiles in terms of a simple sum 
of emission components ignores possible absorption effects. 

According to popular electrodynamic assumptions, the 
accelerating electric field is symmetrical with respect to the 
main meridian. However, the lOF curvature radii are intrin- 
sically symmetric about the caustic phase at a given height 
(and asymmetric with respect to the dipole axis), so cir- 
cular cones centered on this phase are theoretically likely. 
Because of the lOF-curvature asymmetry within the open 
region of fast rotating MSPs, the traditional, dipole-axis- 
centered nest of cones may be considered less natural than 
the decentered geometry of Fig. (5) In the extreme case of the 
fastest known MSPs (P ~ 1.5 ms), the caustic phase occurs 
at the outermost trailing edge of the polar cap, so it be- 
comes difficult to generate a conal structure at all (possibly 
yielding the simplicity of the MP of B1937-I-21). 

The main pulse of J1012-I-5307 exhibits asymmetry of 
flux and width which can be interpreted in terms of the coro- 
tational effects within a radially-extended emission region. 
The observed width of the main pulse suggests that the out- 
ermost components have s ~ 1. This is roughly consistent 
with the widths of the outermost components if r\ ~ i?NS 
and Ar ~ -Rns- However, based solely on the symmetry of 
Piof, the profile should then be symmetrical with respect to 
the caustic phase located somewhere within the main pulse 
window, which is not the case. This implies that either the 
profile is dominated by unknown piof-independent effects, or 
the emission region is located fully on the LS of the caustic 
coordinates. In the uniform-altitude scenario (same r\ and 
r2 for all components) this requires the rings/cones to be 
located asymmetrically with respect to the main meridian 
plane, but the location is not uniquely constrained because 
of the degeneracy of parameters. In the alternative cone- 



altitude scenario, the asymmetry with respect to Scst can 
still be achieved by a region that is symmetrical with re- 
spect to the main meridian in th CF. However, it is hard 
physically to justify the cone-altitude geometry. Neither of 
the scenarios is naturally consistent with the observed width 
of the main pulse, which is the same as the opening angle of 
the near-surface polar beam. 

A possible solution to this problem is a model in which 
the polar tube at low altitudes is fully filled with the ra- 
dio emission, and the minima observed in the MP result 
from absorption effects. The absorption could be provided 
by dense, opaque, radially-extended plasma columns in the 
polar tube. In this model the MP is not regarded as a sum of 
emission components. Rather, it is a single, box-like feature, 
with the minima of different width indented by the plasma 
columns. We expect that the minima will have essentially 
the same widths as the emission components in the emissive 
version of the model presented in Section 2. That is, the min- 
ima on the LS/TS will correspondingly be broad/narrow. 
Then, if the absorption occured (on average) above 0.14_Ric, 
it would be possible to avoid the problem of profile symmetry 
around (j>cBt ■ Such an absorption-based scenario is supported 
by the box-like look of the MP of J1012-I-5307, and presents 
a promising subject for future study. 
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APPENDIX A: EXTENSION OF THE MODEL 
TO NON-ORTHOGONAL PULSARS 

For a 7^ 90° the rotationally-induced acceleration becomes 
an ~ 2ncsinae0 (Dyks 2008) so that the lOF radius of 
curvature becomes 



3 s . 

, ± 2 sm a 



(Al) 



where the plus sign corresponds to the LS, whereas the mi- 
nus is for the TS. 

Near the star surface the equations hold with accuracy 
of ~ 1% (for P ~ 1 s and a as small as 20°). For P ~ 5 ms 
(and a = 20°) the accuracy is degraded down to typically 
~ 10% and can be worse near the locations of zero curvature 
(cf. Fig. |3] for a = 90°). 

If ([^ ~ Q it is also possible to extend the applicability 
of eqs. (O-® to the non- orthogonal case of a 7^ 90°. To do 
this the square root terms in eqs. (HI)-® (but not the linear 
terms) need to be multiplied by 1/ sin a. This is because the 
'not a great circle effect' expands the profile width while 
leaving the linear AR shift intact. 

For example, the detection phase for the TS radiation 
becomes: 

-2-^ + (f,f, (A2) 
2 sin a Y Ric Ric 

so that the caustic amplification determined by the condi- 
tion d^i-v/dr — occurs for radiation emitted from: 

— (A3) 
Ric 64sin2a' ^ ^ 

Thus, for a given r the caustic effects occur more in the 
central parts of the profile (smaller s) than in the orthogonal 
case, whereas for fixed s the caustic enhancement occurs at 
larger r. The amplification is observed at the phase: 

0cst = <^max ^ |r + = 16?1 1 -At- + 0/ , ( A4) 

32 sm a sm a 

ie. the caustic peaks occur further behind the MP than in 
the orthogonal case. 

Though not related to the LS/TS asymmetry, it is also 
easy to find the lOF radius of curvature in the main merid- 
ian plane, where ps -L e^. Using eqs. (O, pil)l. and the 
Pythagorean theorem we get 

/ g 2 \ -1/2 

^-"^HT6;;fc+^^^"'") • 
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